Small world effect in an epidemiological model 
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A model for the spread of an infection is analyzed for 
different population structures. The interactions within 
the population are described by small world networks, 
ranging from ordered lattices to random graphs. For 
the more ordered systems, there is a fluctuating endemic 
state of low infection. At a finite value of the disorder of 
the network, we find a transition to self-sustained oscil- 
lations in the size of the infected subpopulation. 

PACS numbers: 87.23. Ge, 05.65.+b, 87.19.Xx 



I. INTRODUCTION 

How does the dynamics of an infectious disease de- 
pend on the structure of a population? A great amount 
of work has been done on the phenomenological descrip- 
tion of particular epidemic situations 0-0]. A classical 
mathematical approach to these problems deals with well 
mixed populations, where the subpopulations involved 
(typically susceptible, infected and removed) interact in 
proportion to their sizes. With these zero dimensional 
models it has been possible to study, among other epi- 
demic features, the existence of threshold values for the 
spread of an infection || , the asymptotic solution for the 
density of infected people ||-|^], the effect of stochastic 
fluctuations on the modulation of an epidemic situation 
A second classical approach describes spatially ex- 
tended subpopulations. In this, the geographic spread 
of an epidemic can be analyzed as a reaction-diffusion 
process 10 -12|. 

Real populations rarely fall into either of these cate- 
gories, being neither well mixed nor lattices. Recently 
introduced by Watts and Strogatz , small world net- 
works attempt to translate, into an abstract model, the 
complex topology of social interactions. Small worlds 
may play an important role in the study of the influence 
of the network structure upon the dynamics of many so- 
cial processes, such as disease spreading, formation of 
public opinion, distribution of wealth, transmission of 
cultural traits, etc. [O. In relation to epidemiologi- 
cal models, it has been shown that small world net- 
works present a much faster epidemic propagation than 



*E-mail address: abramson@cab.cnea. gov. ar 
^ E-mail address: kupermanQcab. cnea.gov. ar 



reaction-diffusion models or discrete models based on reg- 
ular lattices of a social network |L5| . 

In the original model of small worlds a single param- 
eter p, running from to 1, characterizes the degree of 
disorder of the network, respectively ranging from a reg- 
ular lattice to a completely random graph. It has been 
shown that geometrical properties, as well as certain sta- 
tistical mechanics properties, show a transition at p c = 
in the limit of large systems, N — » oo p6| . That is, any 
finite value of the disorder induces the small world be- 
havior. In this article we show that a sharp transition in 
the behavior of an infection dynamics exists at a finite 
value of p. 



II. EPIDEMIC MODEL 

We analyze a simple model of the spread of an in- 
fectious disease. We want, mainly, to point to the role 
played by the network structure on the temporal dynam- 
ics of the epidemic. The disease has three stages: suscep- 
tible (S), infected (/), and refractory (R). An element of 
the population is described by a single dynamical variable 
adopting one of these three values. Susceptible elements 
can pass to the infected state through contagion by an in- 
fected one. Infected elements pass to the refractory state 
after an infection time rj. Refractory elements return to 
the susceptible state after a recovery time tr. This kind 
of system is usually called SIRS, for the cycle that a 
single element goes over. The contagion is possible only 
during the S phase, and only by an / element. During 
the R phase, the elements are immune and do not infect. 

The interactions between the elements of the popula- 
tion are described by a small world network. The links 
represent the contact between subjects, and infection can 
only proceed through them. As in Watts and Strogatz 
model, the small worlds we study are random networks 
built upon a topological ring with N vertices and coor- 
dination number 2K . Each link connecting a vertex to 
a neighbor in the clockwise sense is then rewired at ran- 
dom, with probability p, to any vertex of the system. 
With probability (1 — p) the original link is preserved. 
Self-connections and multiple connections are prohibited. 
With this procedure, we have a regular lattice at p = 0, 
and progressively random graphs for p > 0. The long 
range links that appear at any p > trigger the small 
world phenomenon. At p = 1 all the links have been 
rewired, and the result is similar to (though not exactly) 
a completely random network. This algorithm should 
be used with caution, since it can produce disconnected 
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graphs. We have used only connected ones for our anal- 
ysis. 

Time proceeds by discrete steps. Each element is char- 
acterized by a time counter Ti{i) = 0, 1, . . . , tj + tr = To, 
describing its phase in the cycle of the disease. The epi- 
demiological state TTi of the element (S, I or R) depends 
on this phase in the following way: 



TTi(t) = S ifr<(t)=0 

7Ti(t) = / if Ti{t) G (l,Tj) 

wi(t) = R if Tj(t) e (rr+1, 



(1) 



The state of an element in the next step depends on its 
current phase in the cycle, and the state of its neighbors 
in the network. The rules of evolution are the following: 

Ti(t+1) = if Ti(t) = and no infection occurs 

Ti(t+1) = 1 if Ti(t) = and i becomes infected 

n(t+i) = n(t)+l if l < n(t) < r 

7i(t + l) = tf n (t)=TR 

(2) 

That is, a susceptible element stays as such, at r = 0, 
until it becomes infected. Once infected, it goes (deter- 
ministically) over a cycle that lasts To time steps. During 
the first tj time steps, it is infected and can potentially 
transmit the disease to a susceptible neighbor. During 
the last tr time steps of the cycle, it remains in state R, 
immune but not contagious. After the cycle is complete, 
it returns to the susceptible state. 

The contagion of a susceptible element by an infected 
one, and the subsequent excitation of the disease cycle 
in the new infected, occur stochastically at a local level. 
Say that the element i is susceptible, and that it has ki 
neighbors, of which fc; n f are infected. Then, i will be- 
come infected with probability fc; n f/&i- Observe that i 
will become infected with probability 1 if all its neigh- 
bors are infected. There are other reasonable choices for 
this mechanism. For example, if the susceptible had a 
probability q of contagion with each infected neighbor, 
we would have a probability of infection [1 — (1 — q) kil11 ]. 
We have tested that both these criteria give qualitatively 
the same results for values of q < 0.2. 



III. NUMERICAL RESULTS 

We have performed extensive numerical simulations of 
the described model. Networks with N = 10 3 to 10 6 ver- 
tices have been explored, with K = 3 to 10. A typical 
realization starts with the generation of the random net- 
work and the initialization of the state of the elements. 
An initial fraction of 0.1 infected, and the rest suscepti- 
ble, was used in all the results shown here. Other initial 
conditions have been explored as well, and no changes 
have been observed in the behavior. After a transient 
a stationary state is achieved, and the computations are 



followed for several thousand time steps to perform sta- 
tistical averages. 

We show in Fig. [I] part of three time series displaying 
the fraction of infected elements in the system, rii n f(t). 
The three curves correspond to systems with different 
values of the disorder parameter: p = 0.01 (top), 0.2 
(middle) and 0.9 (bottom). The three systems have 
N = 10 4 and K = 3, and infection cycles with t/ = 4 and 
Tu = 9. The initial state is random with n;„f(0) = 0.1. 
The 400 time steps shown are representative of the sta- 
tionary state. We can see clearly a transition from an 
endemic situation to an oscillatory one. At p = 0.01 
(top), where the network is nearly a regular lattice, the 
stationary state is a fixed point, with fluctuations. The 
situation corresponds to that of an endemic infection, 
with a low and persistent fraction of infected individuals. 
At high values of p — like the case with p = 0.9 shown in 
the figure (bottom) — large amplitude, self-sustained os- 
cillations develop. The situation is almost periodic, with 
a very well defined period and small fluctuations in am- 
plitude. The period is slightly longer than To, since it 
includes the average time that a susceptible individual 
remains at state 5 1 , before being infected. Epidemio- 
logically, the situation resembles the periodic epidemic 
patterns typical of large populations |3). A mean field 
model of the system, expected to resemble the behavior 
at p = 1, can be easily shown to exhibit these oscillations. 
The transition between both behaviors is apparent — in 
this relatively small system — at the intermediate value 
of disorder p = 0.2, shown in the middle curve. Here a 
low amplitude periodic pattern can be seen, appearing 
and disappearing again in a background of strong fluctu- 
ations. Moreover, the mean value of infection is seen to 
grow with p. 
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FIG. 1. Fraction of infected elements as a function of time. 
Three time series are shown, corresponding to different val- 
ues of the disorder parameter p, as shown in the legends. 
Other parameters are: N = 10 4 , K — 3, n = 4, tr — 9, 
JV inf (0) = 0.1. 

The formation of persistent oscillations corresponds to 
a spontaneous synchronization of a significant fraction of 
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the elements in the system. Their phases r, (f ) in the epi- 
demic cycle become synchronized, and they go over the 
disease process together, becoming ill at the same time, 
and recovering at the same time. We can characterize 
this behavior with a synchronization parameter fj"7| , de- 
fined as: 



a(t) 



1 N 

-Y 



i <pj (t) 



3=1 



(3) 



where c/)j — 27r(rj — 1)/tq is a geometrical phase corre- 
sponding to Tj. We have chosen to let the states r = 
out of the sum in (||), and take into account only the 
deterministic part of the cycles. 

When the system is not synchronized, the phases are 
widely spread in the cycle and the complex numbers e 1 ^ 
are correspondingly spread in the unit circle. In this sit- 
uation a is small. On the other hand, when a significant 
part of the elements are synchronized in the cycle, a is 
large. The synchronization would be strictly a = 1 if all 
the elements were at the same state at the same time. 
However, such a state would end up in Ni n { — after To 
time steps, and the epidemic would end since the system 
is closed and no spontaneous infection of susceptibles is 
being taken into account in the model. 
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FIG. 2. Synchronization of the system as a function of the 
disorder parameter p. Three curves are shown, corresponding 
to different system sizes N = 10 3 , TV = 10 4 , N = 10 5 and 
10 6 , as shown in the legend. Each point corresponds to a 
time average of 2000 time steps, and a subsequent average 
over a number of realizations of the networks and the initial 
condition, as shown in the legend. Other parameters are: 
K = 3, ti = 4, t r = 9, 7V inf (0) = 0.1. 

In Fig. H we show the synchronization parameter cr, 
obtained as a time average of a(t) over 2000 time steps 
and a subsequent average over realizations of the sys- 
tem. Several curves are shown, corresponding to system 
sizes N = 10 3 , 10 4 , 10 5 and 10 6 . A transition in the syn- 
chronization can be seen as p runs from to 1. The 
transition becomes sharper for large systems, at a value 
of the disorder parameter p c sa 0.4. It is worthwhile to 



remark that the transition to synchronization occurs as 
a function of the structure of the network, contrasting 
the phenomenon of synchronization as a function of the 
strength of the interaction, as in other systems of cou- 
pled oscillators Go] . This behavior is observed for a wide 
range of values of 77 and tr. The amplest oscillations 
take place around ti/tr = 1. They disappear when 77 is 
significantly greater or smaller than tr. 

All the systems shown in figures and I have K = 3. 
We have explored higher values of K as well. The picture 
is qualitatively the same, with a sharp transition from a 
quasi-fixed point to a quasi-limit cycle at a finite value 
of p. The critical value p c shifts toward lower values 
for growing K. This is reasonable, since higher values 
of K approaches the system to a globally coupled one 
(at K = N all the elements interact with every other 
one, even at p = 0). So, the mean field behavior (the 
oscillations) can be expected to occur at lower values of 
p. In figure || this effect can be seen for a system with 
N = 10 4 and growing values of K = 3, 5 and 10. 
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FIG. 3. Synchronization of the system as a function of the 
disorder parameter p. Three systems with different values 
of the coordination number K are shown. All of them have 

N = 10 4 , TI =4,TR = 9, JVi„f(0) = 0.1. 



IV. DISCUSSION 

Why does this transition to synchronization take 
place? Unfortunately we do not have analytical argu- 
ments to describe this yet. As mentioned before, a mean 
field model can be shown to predict the oscillations, but 
this can only be expected to describe the network at 
p = 1, and it does not shed light on the nature of the 
transition at lower values of disorder. We can only ad- 
vance some conjectures, based in our observation of the 
dynamical behavior of the system in real time on the 
computer screen. 

An explanation involving L(p), the typical distance be- 
tween pairs of elements (defined as the size of the min- 
imum path connecting two elements) is not plausible, 
since L is known to behave critically with p c = 0, and we 
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observe the transition at p > 0. Complementarity to L, 
small world networks can be described by the degree of 
clusterization C{p) At low values of p, the networks 
are rather regular and highly clustered. As p approaches 
1, C decreases. The crossover from high to low cluster- 
ization occurs at a higher value of p, compared to that 
observed in the decay of L. For systems with K — 3, such 
as those mainly studied here, we have that this crossover 
is mostly concentrated between p — 0.1 and p — 0.5, pre- 
cisely where the onset of oscillations occurs. Moreover, 
the change in the average clusterization C(p) is accom- 
panied by a corresponding one in the distribution of the 
clusterization at the element level, Ci{p). Highly ordered 
networks, at low p, show not only a large value of the 
average clusterization C(p) but also a small dispersion 
around it. On the other extreme, highly disordered net- 
works, with p > 0.5, exhibit a low average clusterization 
and also a low dispersion around it. There is an interme- 
diate range of p, between 0.1 and 0.5, where the average 
clusterization shifts from high to low, but the distribu- 
tions are wide. This indicates that, in this range, the 
system is a mixture of highly ordered, highly clustered 
regions and random, lowly clustered ones. If we consider 
that the clustering structure determines a partition of 
the whole system into smaller, interacting, subsystems, 
the global behavior of the system could be interpreted 
as a superposition of the subsystems' dynamics. When 
p is small the existence of large clusters (essentially one- 
dimensional in our model) inhibits the oscillatory behav- 
ior because, once the infection breaks into such a region, 
it remains restricted to it a long time until all the individ- 
uals have completed the cycle, orderly and deterministi- 
cally. As p increases, the number of big, ordered regions 
decreases. Elements within small regions with some long 
range links go through the infection cycle and become in- 
fected again before long. Some degree of synchronization 
can be seen here (Fig. [y, center, corresponds to p = 0.2). 
When a critical value of p is reached, and the system is 
essentially composed of enough small regions of low clus- 
terization and a similar local dynamics, the synchronized, 
periodic global behavior establishes spontaneously. 

In summary, we have observed a transition at a finite 
value of the disorder in a small world model. The dynam- 
ical behavior of an SIRS epidemiological model changes 
from an irregular, low-amplitude evolution at small p to a 
spontaneous state of wide amplitude oscillations at large 
p. This may be related to observed patterns in real epi- 
demics H , where an effect of the social structure is ob- 
served in the dynamics of the disease. 
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